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Abstract 

The almost Mathieu operator is the discrete Schi'ddinger operator on ^^(Z) defined 

via {Ha,i 3 ,ef){k) = f{k + 1) + f{k — 1) + (3 cos{2Trak +0)f{k). We derive explicit estimates 
for the eigenvalues at the edge of the spectrum of the finite-dimensional almost Mathieu op- 
erator ^ g. We furthermore show that the (properly rescaled) m-th Hermite function pm is 
an approximate eigenvector of ^ g, and that it satisfies fhe same properties fhaf characferize 

' ' (n) 

fhe frue eigenvector associafed fo fhe m-fh largesf eigenvalue of ^ g. Moreover, a properly 

(n) 

franslafed and modulated version of ipm is also an approximate eigenvecfor of ^ g, and if 
satisfies fhe properfies fhaf characferize fhe frue eigenvecfor associafed fo fhe m-fh largesf (in 
modulus) negative eigenvalue. The resulfs hold af fhe edge of fhe specfrum, for any choice of 
9 and under very mild conditions on a and f3. We also give precise esfimafes for fhe size of 
fhe “edge”, and extend some of our resulfs fo Tfa,/ 3 , 0 - The ingredienfs for our proofs comprise 
Taylor expansions, basic lime-frequency analysis, Sfurm sequences, and perlurbalion fheory 
for eigenvalues and eigenvecfors. Numerical simulalions demonsfrafe fhe fighf fif of fhe fheo- 
relical esfimafes. 


1 Introduction 

We consider the almost Mathieu operator on given by 

{Ha, 0 ,ef){x) = f{x + 1) + f{x - 1) + 2/3 cos{ 2 Ttak + 0 )f{x), (1.1) 

with x G Z, /3 G M, a G |), and 9 G [0,27r). This operator is interesting both from a 
phyiscal and a mathematical point of view [fT^[T9ll . In physics, for instance, it serves as a model 
for Bloch electrons in a magnetic field [[T3ll . In mathematics, it appears in connection with graph 
theory and random walks on the Heisenberg group and rotation algebras BUl. A major part 
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of the mathematical fascination of almost Mathieu operators stems from their interesting spectral 
properties, obtained by varying the parameters a, {3, 9, which has led to some deep and beautiful 
mathematics, see e.g [|31|51E1[I51[I71. For example, it is known that the spectrum of the almost 
Mathieu operator is a Cantor set for all irrational a and for all /9 7^ 0, cf. dH. Furthermore, if 
/3 > 1 then exhibits Anderson localization, i.e., the spectrum is pure point with exponentially 
decaying eigenvectors lfT4l . 

A vast amount of literature exists devoted to the study of the bulk spectrum of -ffa,/3,0 and 
its structural characteristics, but very little seems be known about the edge of the spectrum. For 
instance, what is the size of the extreme eigenvalues of how do they depend on a, (3,9 and 

what do the associated eigenvectors look like? These are exactly the questions we will address in 
this paper. 

While in general the localization of the eigenvectors of depends on the choice of (3, 

it turns out that there exist approximate eigenvectors associated with the extreme eigenvalues of 
Ha, 0,6 which are always exponentially localized. Indeed, we will show that for small ot the m-th 
Hermitian function Lpm as well as certain translations and modulations of form almost eigen¬ 
vectors of Ha,0,9 regardless whether a is rational or irrational and as long as the product is 
small. 

There is a natural heuristic explanation why Hermitian functions emerge in connection with al¬ 
most Mathieu operators. Consider the continuous-time version of Ha, 0 ,e in (11.11) by letting x G M 
and set a = 1,/3 = 1, 6* = 0. Then Ha, 0,9 commutes with the Fourier transform on L^(M). It 
is well-known that Hermite functions are eigenfunctions of the Fourier transform, ergo Hermite 
functions are eigenvectors of the aforementioned continuous-time analog of the Mathieu operator. 
Of course, it is no longer true that the discrete Ha,0,9 commutes with the corresponding Fourier 
transform (nor do we want to restrict ourselves to one specific choice of a and 13). But nevertheless 
it may still be true that discretized (and perhaps truncated) Hermite functions are almost eigenvec¬ 
tors for Ha,0,9- We will see that this is indeed the case under some mild conditions, but it only 
holds for the first m Hermite functions where the size of m is either 0(1) or (9(1/y^), where 
7 = iray/^, depending on the desired accuracy (7^ and 7, respectively) of the approximation. We 
will also show a certain symmetry for the eigenvalues of Ha, 0,9 and use this fact to conclude that 
a properly translated and modulated Hermite function is an approximate eigenvector for the m-th 
largest (in modulus) negative eigenvalue. 

The only other work we are aware of that analyzes the eigenvalues of the almost Mathieu 
operator at the edge of the spectrum is ETl . There, the authors analyze a continuous-time model 
to obtain eigenvalue estimates of the discrete-time operator Ha, 0 ,e. They consider the case (3 = 
1, 6* = 0 and small a and arrive at an estimate for the eigenvalues at the right edge of the spectrum 
that is not far from our expression for this particular case (after translating their notation into 
ours and correcting what seems to be a typo in ll2T]| f. But there are several differences to our 
work. First, ll^ does not provide any results about the eigenvectors of Ha,0,9. Second, [|2T]| does 
not derive any error estimates for their approximation, and indeed, an analysis of their approach 
yields that their approximation is only accurate up to order 7 and not 7^. Third, ETll contains no 
quantitative characterization of the size of the edge of the spectrum. On the other hand, the scope 
of ETl is different from ours. 
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The remainder of the paper is organized as follows. In Subseetion ll.il we introduction some 
notation and definitions used throughout the paper. In Section |2] we derive eigenvector and eigen¬ 
value estimates for the finite-dimensional model of the almost Mathieu operator. The ingredients 
for our proof comprise Taylor expansions, basic time-frequency analysis, Sturm sequences, and 
perturbation theory for eigenvalues and eigenvectors. The extension of our main results to the 
infinite dimensional almost Mathieu operator is carried out in Section [3l Finally, in Section |4] we 
complement our theoretical findings with numerical simulations. 


1.1 Definitions and Notation 

We define the unitary operators of translation and modulation, denoted and Mf, respectively by 

(TJ) (x) := fix - a) and (M,/) (x) := e“/(x), 

where the translation is understood in a periodic sense if / is a vector of finite length. It will be 
clear from the context if we are dealing with finite or infinite-dimensional versions of and M^. 
Recall the commutation relations (see e.g. Section 1.2 in (HU) 

TaM, = e-“M6T,. (1.2) 


The discrete and periodic Hermite functions we will be using are derived from the standard 
Hermite functions defined on M (see e.g. (HI) by simple discretization and truncation. We do 
choose a slightly different normalization than in [[T]| by introducing the scaling terms (-y/ 27 )^^ and 
(a/ 27 )^^+^, respectively. 

Definition 1.1 The scaled Hermite functions ipm with parameter > 0 are 


for even m: 


Tmix) 


= 


1=0 


where Cm,i = 


m!(2v/2)2^(-l)f-' 

(20!(f-0! 


(1.3) 


m — 1 
2 


for odd m: pmix) = e 


1=0 


where Cm,i = 


m !( 2 \/ 2 )^^+^(— 1 ) "*2 ^ ^ 

(21 + l)\i^ - l)\ 


(1.4) 

for X G Z and m = 0, 1,.... The discrete, periodic Hermite functions of period n, denoted by 
Tm,n, cire similar to pm, except that the range for x is x = — ..., | — 1 (with periodic boundary 

conditions) and the range for m is m = 0,... ,n — 1. 


(n) 

We denote the finite almost Mathieu operator, acting on sequences of length n, by H)^ ^ q. It 
can be represented by an n x n tridiagonal matrix, which has ones on the two side-diagonals and 
cos(x_^+fc -I- 9) as /c-th entry on its main diagonal for Xk = 2'nak, k = —..., | — 1, and 
j = 0,...,n — 1. Here, we have assumed for simplicity that n is even, the required modification 
for odd n is obvious. Sometimes it is convenient to replace the translation present in the infinite 
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almost Mathieu operator by a periodic translation. In this case we obtain the nxn periodic almost 
Mathieu operator g which is almost a tridiagonal matrix; it is given by 


p("-) 


cos(x_| +9) 1 0 

1 cos(a;_|+i + 6 ) 1 0 

0 1 cos ( a ;_:|_|_2 + 6 ) 1 


0 


1 

0 

0 


0 ... 0 1 COs(x ^_2 + 0 ) 1 

1 0 ... 0 1 cos(xii_i + 6 ) 

' (1.5) 

where Xk = 27rak for /c = — ..., | — 1. If 0 = 0 we write instead of g. 


2 Finite Hermite functions as approximate eigenvectors of the 
finite-dimensional almost Mathieu operator 

In this section we focus on eigenvector and eigenvalue estimates for the finite-dimensional model 
of the almost Mathieu operator. Finite versions of are interesting in their own right. On the 
one hand, numerical simulations are frequently based on truncated versions of on the other 

hand certain problems, such as the study of random walks on the Heisenberg group, are often more 
naturally carried out in the finite setting. 

We first gather some properties of the true eigenvectors of the almost Mathieu operator, col¬ 
lected in the following proposition. 

Proposition 2.1 Consider the finite, non-periodic almost Mathieu operator g. Let > Ci ^ 
• • • > Cn-i be its eigenvalues and ■ ■ ■ ,^n-i the associated eigenvectors. The following 

statements hold: 

= Cmfmfor a// 0 < Trt < 77, - 1; 

2. There exist constants Ci, C 2 , independent ofn such that, for all m 

3. If(r is one of the m-th largest eigenvalues (allowing for multiplicity), then ijjr changes sign 
exactly m times. 

Proof The first statement is simply restating the definition of eigenvalues and eigenvectors. The 
second statement follows from the fact that the inverse of a tri-diagonal (or an almost tri-diagonal 
operator) exhibits exponential off-diagonal decay and the relationship between the spectral projec¬ 
tion and eigenvectors associated to isolated eigenvalues. See (HI for details. The third statement is 
a result of Lemma [Z2] below. ■ 
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Lemma 2.2 Let A be a symmetric tridiagonal n x n matrix with entries ai,... ,an on its main 
diagonal and &i,..., hn-i on its two non-zero side diagonals. Let Xq > X 2 >■■■> \m be the 

eigenvalues of A with multiplicities tq, ..., r^. Let , Vq°\v^\ ..., Vm^'^ be the associated 

ii ) 

eigenvectors. Then, for each 0 < z < m and each 1 < jm < the entries of the vector v] 
change signs m times. That is, for all I < jo < ro, the entries of each of the vectors Vq°^ all have 
the same sign, while for all 1 < ji < ri each of the vectors has only a single index where the 

sign of the entry at that index is different than the one before it, and so on. 

Proof This result follows directly from Theorem 6.1 in |l2l which relates the Sturm sequence to the 
sequence of ratios of eigenvector elements. Using the assumption that bi > 0 for alH = 1,... ,n 
yields the claim. ■ 

The main results of this paper are summarized in the following theorem. In essence we show 
that the Hermite functions (approximately) satisfy all three eigenvector properties listed in Propo¬ 
sition |2T1 The technical details are presented later in this section. 

Theorem 2.3 Let A be either of the operators q or g. Let pm.n be as defined in Defini¬ 
tion o Set 7 = 7iay/fi, let 0 < e < 1, and assume < 7 < ^. Then, for m = 0,..., iV 
where N = 0{1), the following statements hold: 

1. Aprn,n ~ XmTm,n, where Xm = 2/3 + 2e~^ - 4m7e"^; 

2. For each m, there exist constants Ci, C 2 , independent ofn such that 


3. For each 0 < m < n — 1, the entries of change signs exactly m times. 


Proof The second property is an obvious consequence of the definition of Pm,n, while the first and 
third are proved in Theorem 12.61 and Lemma 12.81 respectively. That the theorem applies to both 
^ consequence of Corollary 12.71 ■ 

In particular, the (truncated) Gaussian function is an approximate eigenvector associated with the 
largest eigenvalue of g (this was also proven by Persi Diaconis ifTOll l. 

In fact, via the following symmetry property. Theorem [23] also applies to the m smallest eigen¬ 
values of g and their associated eigenvectors. 


Proposition 2.4 Ifp is an eigenvector of P^^ g with eigenvalue X, then Mi/ 2 Ti/ 2 o‘^ Is an eigen¬ 
vector of P^p g with eigenvalue —A. 


Proof We assume that 9 


P 


(n) 


a,P 


as 


0, the proof for 0 7 ^ 0 is left to the reader. It is convenient to express 

= Ti + T_i + 
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Next we study the eommutation relations between and translation and modulation by eonsid- 
ering 

+ T_i + + (3M_^){TaM,y. (2.1) 

Using (11.21) we have 

and 

Note that = ±1 if & G |Z and = ±1 if aa G |Z. In partieular, if 6 = ^ and a = ^ 

we obtain 

2ck: 2 2ck: 2 

Henee, sinee PjylTj_Miv = —Tj_MiPjylv for any v G C"^, it follows that if ip is an eigenvector 

of P^l with eigenvalue A, then Tj_Mip is an eigenvector of P^l with eigenvalue —A. Finally, 
note that if the entries of p change signs k times, then the entries of Mi 99 change sign exactly 
n — 1 — k times, while the translation operator does not affect the signs of the entries. ■ 

The attentive reader will note that the proposition above also holds for the infinite-dimensional 
almost Mathieu operator Happ. 

In the following lemma we establish an identity about the coefficients in (11.31) and (11.41) . 
and the binomial coefficients 62 Z,fc> which we will need later in the proof of Theorem 12. 6 [ 

Lemma 2.5 For m = 0,..., n — 1; / = 0,..., y, there holds 

‘2Cm,lb2l,2 - 4:Cm,l-lb2l-2,l = -^mCm,l-l, ( 2 . 2 ) 

where bj^k = H) = 

Proof To verify the claim we first note that 62/,2 = and 62 ^- 2 , 1 = 21 — 2. Next, note that 

( 12 . 21 ) is equivalent to 

21 — 2 — m 
^ ^ 2 /( 2 / - 1 ) 

Now, for even m we calculate 

, ( 2 x/ 2 ) 2 T-l)f-^ -8 (f - / + 1 ) m!( 2 \/ 2 ) 2 ^- 2 (-l)f-'+i 

(2 7 = TTl'- _- -_-_ = _—=_ - _ 

■ ( 2 /)!(f-/)! 2 /( 2 /-l)( 2 /- 2 )!(f-/ + l) 

_- 8 (f-/ + l)^ _^2l-2-m^ 

2 /( 2 /- 1 ) ^ 2 /( 2 / - 1 ) 

as desired. The calculation is almost identical for odd m, and is left to the reader. ■ 

While the theorem below is stated for general a, {3, 9 (with some mild conditions on a, (3), it is 
most instructive to first consider the statement for a = 2^/3 = 1,9 = 0. In this case the parameter 
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7 appearing below will take the value 7 = -. The theorem then states that at the right edge of 

(n) (n) 

the speetrum of ipjn,n is an approximate eigenveetor of P^^g with approximate eigenvalue 

Am = 2/3 + 2e~'^ — Am'ye~'^, and a similar result holds for the left edge. The error is of the order 
^ and the edge of the speetrum is of size 0(1). If we allow the approximation error to inerease to 
be of order then the size of the edge of the speetrum will inerease to 0{y/n). 


Theorem 2.6 Let g be defined as in (11.51) and let a, G M+ and 6 G [0, 27r). Set 7 = 
and assume that A < 7 < 1. 

(1) Form = 0,1,..., where N = 0(1), there holds for all x = —..., | — 1 


( Qy9rn,n) (a^) AmV^m,?! ( T 


27ra 


<0(7^ 


(2.4) 


where 
(2) For m 


Am = 2/3 + 2e — 4 m 7 e 

-n + 1, —n + 2,..., —N, where N = 0{1), there holds for all x = — 


(a^) Am( 1) V5m,n T 


9 


27ra 2a 




(2.5) 
- - 1 

( 2 . 6 ) 


where 

Am = -(2,9 + 2e"^ - 4m7e"^). (2.7) 


Proof We prove the result for 6^ = 0 and for even m; the proofs for 9 f) and for odd m are 
similar and left to the reader. Furthermore, for simplieity of notation, throughout the proof we will 
write FI instead of P^’^j and (^m instead of ^rn,n- 

We eompute for x = — ..., | — 1 (reeall that (/?m(f — 1) = </?(—f) due to our assumption of 

periodie boundary eonditions) 

{Hipm){x) = </?m(a; + 1) + Lpm{x - 1) + 2/3 cos(27rQ!x))<yCm(a;) (2.8) 

= e-^"'e-^(e-2T'" ■ 5^Cm,;V(^ + 1)'' + ■ 5^Cm,zV(a^ - 1)'')+ (2-9) 

i i 

+ cos(27rQ!x) (2.10) 

i 


We expand eaeh of the terms (x ± 1)^^ into its binomial series, i.e., 

21 21 

(x+ 1 )^' = and (x- 1 )^' = ^ 

k=0 k=0 


where b 2 i,k = 



( 2 . 11 ) 


7 









and obtain after some simple ealeulations 


m/2 


/=0 


= e ^b 2 i,oie + 2/3cos(27raa;) 


ml2 


21 


1=0 k=l 

= (I) + (11). 


( 2 . 12 ) 


(2.13) 


(2.14) 


We will now show that (I) = (pm{x) ■ (2/3 + 2e“'^) + 0(7^) and (II) = —ipm(x) ■ 4 m 7 e“^ + C>( 7 ^), 

from whieh (12.41) and (12.51) will follow. 

We first eonsider the term (I). We rewrite (I) as 


(I) = e- 


— e-'y^ I 


+ + 2/3 cos(27raa:) j ^ 


m/2 


Cm,n^x‘^K 


1=0 


Using Taylor approximations for , and cos(27rax) respeetively, we obtain after some 

rearrangements (whieh are justified due to the absolute summability of eaeh of the involved infinite 
series) 




('g- 7 (g- 27 X ^ 2/3 cos( 27 rax)) = 


3! 


(2.15) 




3! 


= ( 2 e-^ + 2 e-^ 


- 7 ( 273 ;)' 


(27rQ;x)^ 


21 + 2/3 — 2/3--he (-Ri(a;) + i?2(3;)) + 2/3i?3(a:)j, (2.16) 

where Ri , R 2 , and Ro are the remainder terms of the Taylor expansion for and cos(27rQ!a;) 

(in this order) respeetively, given by 


R^ (x) = 


(-27)^e-2T'«i 


X'^, R2(x) = 


(27)^e2T'«2 


x^, Rsix) = 


(2™)"^ cos(^3) 


X 


4! 1 , ov 7 

with real numbers ^ 1 , (^ 2 , ■Cs between 0 and x. We use a seeond-order Taylor approximation for e“'*' 


with eorresponding remainder term i? 4 ( 7 ) = ^ 7 ^ 7 ^ (for some .^4 G ( 0 , 7 )) in (12.161) . 


Henee (12.161) beeomes 


e ( 2e + ( 270 :)^ + ( -7 + ^ + Ri{'y) ) ( 270 ;)^ + 2/3 - f3{27raxy ) + 


+e 


—^X 


1 — 7 + (^ + i? 4 ( 7 )^ (.^ 1 ( 3 ;) + R 2 {x)) + 2(3R^{x)^ . 


(2.17) 

(2.18) 
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Since o' = na^/B. the terms (2'yx)^ and —B(27iaxV in (12.17b caneel. Clearlv. |i?.i (a;) I < 

|i? 2 (a^)| < and 1 7 ^ 4 ( 7 ) I < It is convenient to substitute a = in Rz{x), in which 


4 ! 


case we get |i? 3 (a:) | < . Thus, we can bound the expression in (12.18b from above by 


„2 


(( - 7 + |- + ^ 4 ( 7 )) (fli(i) + R^S^)) + 2m3(.x)) 


<^-’'((7 + ^ + ^)( 


7 ^ ( 7 )^, /( 270 :)“^ e^'^^{2'yx) 


4! 


+ 


4! 


+ 2/3 


( 273 ^) 

/^M! 




(2.19) 

( 2 . 20 ) 


Assume now that |x| < ^ then we can further bound the expression in (12.20b from above by 


^-■yx- 


7 ^ 7^4 4 ( 270 :)'^ e‘^'^^{2'yxy 


«'W) 


4! ^ 4! PH\ 


(h+ir+ i!)( 

Moreover, if |x| < we can bound the term ( — 7 + ^ + 774 ( 7 )) in (12.17b by 


( 2 . 21 ) 

( 2 . 22 ) 


(-77| + fi4(7))^^ 


<272 +73 + 1. <0(^2). 


Now suppose |x| > -. We set |x| = ( 7 /- ) for some c with 


1 < c < 


2 log(n/ 2 ) 
log(l/ 7 ) ' 


(2.23) 


The upper bound in (12.23b ensures that |a;| <7 and the assumption ^ < 7 < 1 implies that 
condition (12.23b is not empty. Then 




I ( 7 ) 3 , 4 ( 270 :)'^ e^'^^{ 2 'yx) 


“ (h + i!+^)( 4! 


, ( 27 X) 


+ + - 

4-2c 2y^-<^ fy„,\A-2c a^A-2c 


'(27) 


4-2c 


+ ^^))<0{A), (2.25) 


(2.24) 


4! 


3/3 


where the last inequality follows from basic inequalities like e~'^^ ^ ^2 

follows from for all y >0). 

Thus we have shown that 

e- 7 + j'g-7(e-27^ + e^^^) + 2/3 cos(27raa:)) = (2e-3' + 2/3) + (2.26) 

Returning to the term (I) in (12.12b . we obtain, using (12.26b . 


m/2 


(I) = (/5m(a:)(2e ^ + 2/3) + ((7.^ ^ Cm,n^x'^^ 


(2.27) 


i=0 
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where = C>( 7 ^). 

Let us analyze the error term Y1T=q Cm,n^x^K Using Stirling’s Formula ( m Page 257]) we 
note that Cm,i grows at least as fast as but not faster than (m/e)'^. Henee, the error term 

will remain of size £>( 7 ^), as long as we ensure that m does not exeeed C>( 1 ). 

We now proceed to showing that (II) = — 99 ^( 2 ^) • (4m7e“'^) + 0 ( 7 ^). The key to this part of 
the proof is the observation that H acts “locally” on the powers that appear in the definition of 
Lfm- Recall that the term (II) has the form 


m/2 21 

1=0 k=l 


(2.28) 


Analogous to the calculations leading up to (12.261) we can show that for odd k there holds 

e-7-^e-7(e-27x + (2.29) 

and for even k 

e-7^^e-^(e-2'^^ + (-1)^6^^^) = + 0 ( 7 ^). (2.30) 

Using (12.291) and (12.301) we can express (12.281) as 


(g-7^^g-7 


ml2 


(2.31) 


^i.l‘^))\^(^m,ll\'^{,-^lx)h2l,kX^^ ^ 2,1)21,kX^^ 

/=0 odd k even k 

Furthermore, using estimates similar to the ones used in deriving the bounds for (I), one easily 
verifies that 

g- 7 a;^g- 77 ^(_ 472 ;) 3 ; 2 i-fc ^ 0 ( 7 ^), for /c > 3, 


and 

Therefore, 


< 0 ( 7 "), 

mf2 


for /c > 4. 


g-7^ g-7 


^Cm,n^{^{-^'yx)b2l,kx'^^ ^2b2l,kx‘^^ 


1=0 


odd k 

m/2 


even k 


(e + 0 ( 7 ^)) ^c^,z 7 '((- 47 a:) 62 z,ia;^' ^ + 262 ^, 22 ;^' 


1=0 


The expression above implies that H acts “locally” on the powers of Lfm- 
Moreover, 


(2.32) 

(2.33) 


m/2 


g-7^'g-7j 

y/ c^,z7/(-473^)^2Z,i^^^ ^ + 262z,2a;^' 

/=0 

(2.34) 

g-7^^g-7 

“(-47m)c^,^7"^/2^"* + 2cm,f7"^^^bm,2X^-^ 

(2.35) 

m/2—] 




{-i'j)cm,n'’b2i,ix‘^'’ + 2 Cm,z 7 ' 62 Z, 2 X^'"^) , 

(2.36) 


^=0 
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and 


m/2—1 


^ ^ {-4:'y)cm,n^b2i,ix^^ + 2cm,n^h2i,2X^^ 


1=0 


ml2 


(-4m7)cm_-7”"/^x™ + ^ {2cm,a^b2i,2X^’' ^ - 47Cm,;_i7^ %i-2,ix‘^^ 


i=0 


Tnl2 


(-4m7)cm,-7”"/^x™ + 7 ^(2cm,/fc2Z,2 - 4cm,z-i&2z-2,i) 


z=o 


m/2 


= (-4m7)cm,f - 4m7 y/ c^,z-i7 

m/2 


1-1^21-2 


1=0 


4m'yy^Cm,n^x‘^\ 


(2.37) 

(2.38) 

(2.39) 

(2.40) 


i=0 


where we have used Lemma [231 in (12.391) . Henee, up to an error of size (9(7^) we have 


g-7^ £-7 


m/2 


m/2 




2Z,fca; 


2Z-Zc 


y/ ‘^b2i,kx‘^^ = -4m7e "^e ^ Cm,z7 




«=o 


odd k 


even k 


1=0 


(2.41) 

Analogous to the estimate of (I), we need to ehoose m to be not larger than (7(1), in order to 
keep the error term 


m/2 


0{'y‘^)^'^Cm,n^{{-^'yx)b2i,ix'^^ ^ + 2b2l,2x‘^^ 
1=0 


(2.42) 


of size (7(7^). 

Therefore, for sueh an m, by invoking (11.31) . (12.311) . and (12.411) we can express the term (II) in 
(12.281) as 

(II) = -4m7e"'^v5^(a:) + (9(7^). 

Hence we have shown that 

{Hipm){x) = (2/3 + 2e"'^ - 4:mje~'^)(pm{x) + (7(7^), 

which establishes claims (12.41) and (12.51) . 

Claims (12.61) and (12.71) of Theorem 12.61 follow now from Proposition [2]4l ■ 


Remark 2.1 If we allow the error bound in (12.41) and (12.61) to increase from (7(7^) to ( 7 ( 7 ), then 
the range for which the ifm serve as approximate eigenfuctions for g extends from m (7(1) 

to m ^ 0{^JYp/). This can be seen by a suitable modification of the proof of Theorem 12.61 m the 
derivations in (12.241) . (12.251) and in the error term (12.271) . as well as in the calculations leading up 
to (12.411) and in the error term (12.421) . (We leave the details to the reader.) See also the numerical 
simulations depicted in Figure^ which illustrate this fact. 
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Corollary 2.7 Theorem \2.6\ also holds if we replace ^ with g. 


Proof A simple application of the triangle inequality yields 


T-jiX^') \ 


< 


^a,0,e^rn,n ^aPfi^rn^r 


+ 


TD (^) \ 


We compute 


^ap,e^rn,n ^a,l3,0^m,n 




n 


- 1 


m/2 


i=0 


+ 


+ 

m/2 


'pr. 


n 


^ :7vq-? 


E' 

z=o 


21 


(2.43) 

(2.44) 


Now, from 12.27] we know that the sums are of size C>( 7 ^), so we conclude (in a large overestimate 
when n is large) that 

— PaPfiPm,n ^ = 0{') ). 

Meanwhile, from Theorem 12.61 we know that 




Lemma 2.8 Let pm,n be as before. Let 0 < e < 1 and assume < 7 < ^- Then, as long as 
m < rf — 1, the entries of prn,n change signs m times. 

Proof Recall that the Hermite polynomial of order m has m distinct real roots. From |]T^ . we 
know that the roots of the mth Hermite polynomial lie in the interval [—\/2m + 2, ^/2mPr2]. 
Now, in our definition of Prn,n^ we have applied the transformation x i—)■ y/^x to the Hermite 
polynomials. Thus, in our case, we see that the zeroes of prn,n lie in the interval 


7 


m + 1 jm + 1 


7 


Now, the condition < 7 implies that the zeroes of prn,n lie in the interval 


1_£ 1_£ 
n 2 - n 2 . - 

Vm + 1, Vm + 1 


from which we see that so long as m < — 1, all of the zeroes of the mth Hermite polynomial 

lie in the interval 

■ n n' 

.~ 2 ’ 2 . ■ 
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Then, if the distance between consecutive zeroes is always larger than 1, the lemma will follow. 
Now, from Il20l we know that the distance between consecutive roots of the mth Hermite polyno¬ 
mial is at least . Applying our scaling, we see that the distance between consecutive zeroes 
of ifm n is at least , ^ The condition 7 < ^ leads to 

’ y/2'f{2m+l) "■ 

71 I rA 

^27(2m + 1) ^ 4m + 2’ 

and so we see that whenever m < , the minimum distance between consecutive zeroes of 

^m,n will be greater than 1. Finally, note that — 1 < 

■ 

In particular, by choosing e = 4 in Lemma [2^ and thus < 7 < we obtain that for 
m < ^/n + 1, (pm,n changes sign m times. 

Remark 2.2 Let a = ^, where r G N+. Then it is easy to see that each eigenvalue of has 

multiplicity (at least) ^ and if f is an approximate eigenvector of with eigenvalue A, then so 

isTkn/rfJork = 1 ,... ,r - 1 . 

3 Eigenanalysis of the infinite almost Mathieu operator 

In the preceding sections we have seen that prn,n are approximate eigenfunctions with associated 
eigenvalues Am for both P^^q and We now consider the infinite dimensional operator 

and its relationship with the non-periodic version of the Hermite functions defined in (11.31) 
and (11.41) . We shall denote these functions by pm, and for all m > 0 they are defined exactly as 
in Definition 11.11 except that the range of x is Z. Note that the claims in Proposition 12.li no longer 
make sense for the infinite dimensional operator Hapfi- Nevertheless, we see that Am and pm still 
function as approximate eigenpairs for Hap^, and the statements of Theorem 12.31 still hold. 

Theorem 3.1 Let Happ, Pm, Am be as before. Set 7 = Tia^/P and assume 0 < 7 < 1. Then, for 
m = 0,..., A^, where N = 0(1), the following statements hold: 

1. Hq^ppPjyi ~ SraPm,) 

2. For each m, there exist constants Ci,C 2 > H such that 

\Tm('i)\ < 

3. For all 0 < m < ^ — ^, the entries of pm change signs exactly m times. 

Proof The first statement is proved in Corollary 13.21 to Theorem 12.61 below. The second follows 
from the definition of the Hermite functions. The third statement is an easy consequence of Lemma 

m m 
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Corollary 3.2 Let be as defined in (11.11) and let a, /9 G M"*" and 6 G [0, 2ti). Set 7 = 
and assume that 0 < 7 < 1 . 

(1) Form = 0,1,..., A^, where N = 0{1), T-e ipm is an approximate eigenvector for Ha,p ,9 with 

277 0: 

approximate eigenvalue = 2/3 + 2e~'^ — Moreover, there holds for all x E Z 


{,bia,P,dFrn){pb) Xm T-6 Prn.(x) 

2770 




(3.1) 


(2) For m = —n + 1, —n + 2,..., —N, where N = 0(1), T-e is an approximate eigen- 

2770 

vector and \m = —(2/3 + 2e“'>' — is an approximate eigenvalue associated with the m-th 

largest negative eigenvalue of Ha,p- Moreover, there holds for all x eZ 


(Ha,P,dFm)(^) I’^rnT n — e Mil <~Prn. (x) 

2770 2 


< o(^^). 


(3.2) 


If we replace 0(7^) in (13.11) and (13.21) by 0('y), then the range for m increases from N = 0(1) to 

/v = 0(717). 


Proof As in the proof of Theorem I2.6[ we prove the result for even m and 9 = 0. Indeed, the 
proof proceeds exactly as before, except that in this case we no longer need the assumption (12.231) . 
Thus, (13.11) holds for 0 < 7 < 1 and for all x E Z. We apply the infinite-dimensional version of 
Proposition |2]4] to obtain (13.21) . ■ 


4 Numerical simulations 

(n) 

In this section we illustrate our theoretical findings with numerical examples. We consider Pa 
with the parameters n = 10000 , a = ^, (3 = 1,9 = 0, hence in this case 7 = 7 . Figure [T| shows 
the true eigenvalues of and the approximate eigenvalues given by formula (12.51) . Depicted are 
the 600 algebraically largest eigenvalues and their approximations. 

Theorem [2]6] predicts that the approximation error will be of size 7 ^ for the first few approxi¬ 
mate eigenvalues, while about y/n many approximate eigenvalues deviate from the true eigenvalues 
by an error of less than about 1/n. Figure |2]confirms this prediction. In Figure |2] we show the num- 
ber of approximate eigenvalues that approximate the true eigenvalues of ^ ^ within an error of 7 
(red solid line) as well as within an error of 7 ^ (green dash-dotted line). Also depicted is the graph 
of 1/ y/y (blue dashed line). The picture confirms clearly our theoretical findings, which imply that 
the number of approximate eigenvalues that deviate from the true eigenvalues by an error less than 
7 ^ is of size 0 ( 1 ), i.e, independent of n, while the number of approximate eigenvalues that deviate 
from the true eigenvalues by an error less than 7 scales like 1 / y/y, thus in this example like 
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True and approximate eigenvalues 




Figure 2: Shown are the number of approximate eigenvalues that approximate the true eigen¬ 
values of Q within an error of 7 (red solid line) as well as within an error of 7^ (green 
dash-dotted line). Also depicted is the graph of 1 /^/7 (blue dashed line). The values for ^ 
are n = 1000,... , 30000, a = ^, /3 = 1,0 = 0, hence in this example 7 = 
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